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Abstract 


Calculations are presented for electric quadrupole excitations in relativistic nucleus- 
nucleus collisions. The theoretical results are compared to an extensive data set and it is 
found that electric quadrupole effects provide substantial corrections to cross sections, 
especially for heavier nuclei. 


L INTRODUCTION 


The search for a fundamentally new state of matter in the form of a Quark-Gluon 
Plasma *) has stimulated the production of very high energy nuclear beams. The hope is to 
observe the Quark-Gluon Plasma in a relativistic nucleus-nucleus collision. At the 
Berkeley Bevalac a variety of light nuclei such as 12 C, 1 ^0 and 20 Ne can be accelerated up 
to energies of 2. 1 GeV/N and heavier nuclei such as 139 La and 238 U can be accelerated to 
1.26 and 0.96 GeV/N respectively. At Brookhaven, New York, 16 0 beams are available at 
14.6 GeV/N and at the CERN SPS, in Geneva, beams of 1 ^o and 32 S are both produced at 
60 and 200 GeV/N. The Relativistic Heavy Ion Collider (RHIC) is expected to produce 
two colliding beams at 100 GeV/N to give a total center-of-mass energy of 200 GeV/N, 
which corresponds to a single beam energy of 21 TeV/N. Grabiak 2 ) has pointed out that 
nuclear beams of 3.5 TeV/N and 8 TeV/N may be possible at the CERN Large Hadron 
Collider (LHC) or the Superconducting Super Collider (SSC). By way of comparison, the 
majority of Galactic Cosmic Rays have energies 3 ) of about 1 GeV/N, with a range 
typically from 10 MeV/N to 1 TeV/N. However, the JACEE collaboration 4 ) has made 
observations as high as 1000 TeV/N. 

Nucleus-nucleus reactions proceed mainly through either the Strong or Electromagnetic 
(EM) interactions. Historically, Strong interaction processes have been the main object of 
study 5 ), however with the availability of the above high energy nuclear beams there has 
been a resurgence of interest in EM interactions in relativistic nucleus-nucleus collisions. 6) 

The primary theoretical tool for studying these relativistic EM processes has been via 
the Weizsacker-Williams (WW) method 6 ‘ 7 ) of virtual quanta. The nucleus-nucleus total 
EM reaction cross section is 

Ct = j Nww (Ey) tf(Ey) dEy (1) 

where Ey is the virtual photon energy, N^y (Ey) is the WW virtual photon spectrum and 
c(Ey) is the photonuclear reaction cross section. For high accuracy it is important to use 
experimental photonuclear data for o(Ey). (For an excellent compilation of photoneutron 


data see reference 8.) However, a more exact formulation of a involves a breakdown into 
the various EM multipolarities such as electric dipole (El), electric quadrupole (E2), 
magnetic dipole (Ml) etc. The most important contributions to a are from El and E2 so 

that 


a - a m + a E2 


I 


[NEi(Ey)OEi(Ey) + N E2 (E y )a E 2(E y )] dEy 


( 2 ) 


where Ngi (Ey) is the virtual photon spectrum of a particular multipolarity due to the 
projectile nucleus and Ogi (Ey) is the photoriuclear reaction cross section of the target 
nucleus. Bertulani and Baur 6) have derived expressions for Ngi (Ey) and found that the 
electric dipole spectrum is the same as the WW spectrum, i.e. Ngi (Ey) = N w (Ey). 
Furthermore at very high projectile energies all Ngi (Ey) and Nmi (Ey) are equal so that 
equation (1) is seen to be a very high energy approximation to all multipolarities included in 
equation (2). Bertulani and Baur 6 ) have made a crude estimate of the EM cross section 
using equation (2) but they pulled Nei (Ey) and Ne 2 (Ey) outside the integral and evaluated 

them at a single energy and used sum rules to evaluate J °Ei (Ey) dEy. However, this 

procedure did not yield very accurate results. Thus I undertook a more exact study 9 ) 
leaving equation (2) as it stands, and using experimental data for the photonuclear cross 
sections by defining 

(3) 

^E1 (Ey) = <3 eX pt. (Ey) — <TE2(Ey) 

where a ex pt. (Ey) is the experimentally measured photonuclear cross section and (Te 2 (Ey) is 
a theoretical calculation based on a Lorentzian shape for the Electric Giant Quadrupole 
Resonance (GQR). Details for this procedure can be found in reference 9. As was noted 
in that reference, the above procedure yields very accurate values for the sum agi + cj£2 
(which is to be compared to nucleus-nucleus reaction experiments) even though the GQR 
parameters are uncertain. The basic reason for this, as can be seen from equation (3), is 


that an under (over) estimate in Oe 2 (Ey) will give an over (under) estimate in agi (Ey), so 
that the combined <7 ei + <5E2 in equation (2) will not change very much. 

In reference 9 a detailed study of El and E2 was undertaken for the reaction 89 Y (RHI, 
X) 88 Y where RHI refers to various Relativistic Heavy Ions and X is anything. It was 
found that E2 effects account for a considerable fraction of the cross section, and that 
inclusion of E2 (via equation 2) provides improved agreement with experiment over the 
WW method. Given this situation, it was decided to compare this theoretical approach to 
as much experimental data as possible. Thus the present work involves a comparison to 
neutron emission from 89 Y, 197 Au, 59 Co and neutron and proton emission from 12 C, 16 0 
and 18 0 which includes both electric dipole and quadrupole effects. This complements 
earlier work 7 ) which involved an extensive comparison of WW theory to experiment. 
CALCULATIQNAL METHOD 

The basic calculational method is outlined in reference 9 and the discussion will not be 
repeated here. Also reference 7 includes a very detailed summary of which photonuclear 
data were used for a expt . (Ey) in equation (3). The same data is used in the present work. 
All isoscalar GQR parameters were taken from the compilation of Bertrand 10 ) and are 
listed in Table 1. As mentioned in the Introduction, even though these parameters are 
somewhat uncertain the total EM cross section ogi + 0£2 is expected to be very accurate 9 ) 
due to the subtraction procedure of equation (3). The most inaccurate results would be 
expected for the 12 C, 16 0, 18 0 GQR parameters where the isoscalar GQR is fragmented 
into several components 10). Only a single Lorentzian 9 ) was 



Table 1: Isoscalar Giant Quadrupole Resonance (GQR) Parameters taken from the 
compilation of Bertrand 10 \ E is the GQR resonance excitation energy, T is the full-width 
at half maximum and f is the fractional depletion of the Energy Weighted Sum Rule. (The 
GQR of light nuclei are fragmented into several peaks, so that the parameters below 
represent an estimated average value.) 


Nucleus 

E 

(MeV) 

r 

(MeV) 

f 

12C 

27. 5 a 

3.0 b 

0.4 d 

16o 

22.0 b 

3.0 C 

0.4 d 

18o 

24.0 3 

3.0 C 

0.4c 

59Co 

16.3 b 

5.6 b 

0.61 b 

89y 

13. 8 b 

3.2 b 

0.55 d 

192 Au 

10.8 a 

2.9 b 

0.95 d 


a E is calculated from 63 A' 1 / 3 
b best value taken from Table 4 of ref. 10 
c estimate 

d from Fig 23 of ref. 10 

used in the present work. However oe 2 is found to be quite small for these nuclei (see 
below) so that my conclusion that the calculated cjei + OE2 is accurate remains valid. 

For the nuclei 12 C, 16 0 and 18 0, proton (p) emission occurs as well as neutron (n) 
emission. Thus equation (3) needs to be modified to incorporate the branching ratio. I 
assume that the excited nucleus decays only by proton or neutron emission and that the 
(photon) energy dependent neutron branching ratio is defined as 

^expL(Ey>n) + ^expLCEyiP) 


so that 


(5) 


CE2 (Ey, n) = f n (Ey) OE2 (Ey) 

where Ge2 (Ey) is the photonuclear GQR cross section. Thus for proton and neutron 
emission equation (3) becomes 

Cei (Ey, n) = Gexpt. (Ey, n) - f n (Ey) ge2 (Ey) (6a) 

and 

Oei (Ey, p) = G exp t. (Ey, p) - [1 - f n (Ey)] Ge2 (Ey) (6b) 

Equations (4) - (6) were used for nucleon emission from 12 C, l6 0 and 18 0. For 59 Co, the 
(y, p) cross section is not available and so a constant value of f n = 0.7 (suggested from 
reference 11) was used. For 89 Y and 197 Au I used f n = 1.0. 

RESULTS AND DISCUSSION 

The calculated results are listed in Table 2, along with the experimental results of 
various groups. 12-16) <j E1 + <j E2 is the calculated result to be compared with the data 
Gexpt- Also listed are the results of WW calculations. 7 ) In all cases two theoretical cross 
sections are listed. The first is calculated using an expression for the minimum impact 
parameter as 

b m in = Ro i(T) + R 0 i (P) (7) 

where Roi represents the 10-percent charge density radius 7 ) of the target or projectile. 

The second theoretical cross section listed in parentheses in Table 2 uses b m j n given by Hill 
et al. 14-16) ^ 

bmin = r 0 [A p V3 + A T W - X(A p l/3 + Ar 1 ' 3 ] (8) 

where r 0 = 1.34 fm and X = 0.75. (Note that my WW calculations disagree with earlier 
results of Hill et al. 14-16) due to an error in their calculations. 19-20)) 

There are several features readily apparent from Table 2; 

(0 Gei + Ge 2 is always larger than Gww- However, for nucleon emission from 12 C, 
16 0 and 18 0 this difference is never larger than about 4%, but for neutron emission from 
59 Co, 89 Y and 197 Au the difference is much larger varying between about 7% - 15%. 



“) For nucleon emission from 12 C and 16 0 both oei + <*E2 and agree with 
experiment for both choices of bm^. 

“0 For nucleon emission from 18 0 both ctei + GE2 and disagree with experiment 
for both choices of bmin. y actually gives slightly better agreement but not by a 
significant amount. 

iv ) For neutron emission from 197 Au, Gei + <*E2 is significantly closer to experimental 
values than is Oww> although for most cases it still lies outside the error bars. An 
exception however is a much poorer agreement for 139 La (see also references 19 and 20). 
Significant discrepancies with 197 Au data have been noted previously for WW theory. 7 ) 

v ) For neutron emission from 89 Y, oei + &E2 is in much better agreement with 
experiment than is oww- This is especially true for the 40 Ar and 56 Fe projectiles. 

vi ) For 59 Co, oei + <TE2 is again better for 20 Ne, although slighdy worse for ^Fe. As 
above the agreement for the 139 La projectile is significantly poorer. 

S U MMARY AND CO NC LUSIO NS 

Calculations have been made for nucleon emission via EM dissociation in relativistic 
nucleus-nucleus collisions. Results are presented for Weizsacker-Williams theory and also 
for separate electric dipole and quadrupole components. The theories have been compared 
to an extensive data set. It is found that electric quadrupole (E2) effects are not significant 
for proton and neutron emission from 12 C, 16 0 or 18 0. However, E2 contributions are 
substantial for neutron emission from 39 Co, 89 Y and 197 Au, generally leading to improved 
agreement between theory and experiment. Notable disagreements occur for 139 La 
projectiles (1.26 GeV/N) where the theoretical oei + &E2 are too big. Quadrupole effects 
improve the theoretical results for 16 0 projectiles at 60 and 200 GeV/N, although the 
theoretical cross sections are still too small. 

In general it has been found that electric quadrupole effects are an important component 
in nucleus-nucleus collisions and that these effects can be calculated accurately. 
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Table 2 Calculated results, cei + Oe 2 and cww> compared to experiment. 12-16) t wo theoretical cross sections are listed, 
first set uses b m i n given by equation (7) and the second set (in parentheses) uses b m i„ given by equation (8). 



15n I 30±16 | 42 (46) | 42 (46) | 1 (2) | 43 (48) 




" n B -1±4 0.6 (0.7) 0.6 (0.7) 0(0) 0.6 (0.7) 

1.05 U C -2±5 0.3 (0.4) 0.3 (0.4) 0(0) 0.3 (0.4) 

" n B -2±5 0.5 (0.6) 0.5 (0.6) 0(0) 0.5 (0.6) 
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Table 2 (cont.) 
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